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FOREWORD 


This  report  describes  a  portion  of  the  results  obtained  on 
NASA  Grant  NSG  3044 .  This  work  was  done  under  subcontract  to  the 
University  of  Illinois,  Urbana,  with  Prof.  S.S.  Wang  as  the  Principal 
Investigator.  The  prime  grantee  was  the  Massachusetts  Institute  of 
Technology,  with  Prof.  F.J.  McGarry  as  the  Principal  Investigator 
and  Dr.  J.F.  Mandell  as  a  major  participant.  The  NASA  -  LeRC  Project 
Manager  was  Dr.  C.C.  Chamis . 

Efforts  in  this  project  are  primarily  directed  towards  the  de¬ 
velopment  for  finite  element  analyses  for  the  study  of  flaw  growth  and 
fracture  of  fiber  composites.  This  report  presents  exact  solutions  for 
edge  delaminations  in  angle-ply  composites. 
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ABSTRACT 


Edge  delamination  has  caused  severe  concern  in  the  design  and 
analysis  of  advanced  composite  materials  and  structures.  Due  to  its 
complex  nature,  very  limited  knowledge  for  the  problem  is  currently 
available.  It  involves  not  only  geometric  and  material  discontinuities 
but  also  inherently  coupled  mode  I,  II  and  III  fracture  in  the  layered 
anisotropic  system.  Based  on  complex-variable  stress  potentials  in 
the  anisotropic  elasticity  theory  and  eigenfunction  expansion,  exact 
orders  of  the  crack-tip  stress  singularity  and  complete  field  solutions 
are  obtained.  Results  are  given  for  edge-delaminated  composites 
subjected  to  uniform  axial  extension  for  illustrative  purposes. 

Effects  of  geometric,  lamination,  and  crack  variables  are  determined. 


1 .  INTRODUCTION 


Edge  delamination  is  frequently  encountered  in  angle-ply  compo¬ 
site  laminates.  It  is  due  to  high  stress  concentrations  at  geometric 
boundaries  and  the  inherently  weak  interlaminar  strength  along  the  ply 
Interface.  The  problem  posed  by  edge  delamination  is  of  great  theoreti¬ 
cal  interest.  It  also  is  of  significant  technical  Importance  in  determining 
the  structural  integrity  and  damage  tolerance  of  advanced  fiber-rein¬ 
forced  composites  and  their  applications  to  advanced  engineering 
structures  and  components.  The  presence  and  growth  of  delamination 
cracks  from  geometric  boundaries  of  composite  laminates  may  lead  to 
severe  reliability  and  safety  problems  of  fiber  composite  materials 
and  structures  such  as  the  reduction  of  structural  stiffness,  the 
exposure  of  the  interior  to  adverse  environmental  attack,  and  the 
disintegration  of  the  material,  which  may  cause  the  final  failure. 

Thus,  understanding  the  basic  nature  of  edge  delamination  is  of 
critical  importance  in  damage  characterization  and  accurate  assessment 
of  flaw  criticality  and  structural  integrity  of  advanced  composites. 

The  edge  delamination  problem  is  very  complex  in  nature  and 
extremely  difficult  to  solve  •  It  involves  geometric  and  materials 
discontinuities,  i.e.,  free  edges,  interlaminar  cracks  and  variation  of 
ply  properties  in  the  transverse  direction.  It  also  involves  inher¬ 
ently  coupled  mode  I,  II  and  III  fracture  in  the  anisotropic  layered 
material  system  such  as  angle-ply  composite  laminates.  Edge  de- 
lamlnatlon  is  basically  a  fracture  problem  involving  an  interfacial 
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crack  between  two  highly  anisotropic  fiber-composite  laminae 
under  general  loading  conditions.  The  problem  of  an  Interfacial  crack 
between  two  dissimilar  isotropic  materials  has  received  much  atten¬ 
tion  recently,  for  example.  Refs.  [1-12].  But  the  study  of  a  de¬ 
lamination  between  strongly  anisotropic  fiber  composite  layers, 
especially  in  finite-dimensional  laminates  under  general  loading, 
has  been  very  limited,  to  the  author's  knowledge.  In  this  paper,  a 
rigorous  investigation  of  the  basic  nature  of  the  coupled  opening, 
sliding  and  antiplane  shearing  fracture  behavior  of  edge  delamination 
is  presented  for  composite  laminates  under  uniform  axial  extension. 
Basic  formulation  of  the  problem  based  on  the  theory  of  anisotropic 
elasticity  and  eigenfunction  expansion  is  given  in  the  next  section. 
General  solutions  and  associated  stress  singularities  for  the  edge 
delamination  are  derived  in  Section  3.  Fracture  parameters  such  as 
mixed-mode  stress  intensity  factors  K^.,  and  as  well  as  the 

energy  release  rate  G  are  defined  and  determined  in  Section  A.  Num¬ 
erical  results  for  edge  delaminated  composite  laminates  subjected 
to  uniform  axial  loading  are  shown  in  Section  5  to  illustrate  the 
fundamental  behavior  of  edge  delamination  cracks.  Effects  of  geo¬ 
metric,  lamination  and  crack  variables  are  studied  in  detail. 
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2 .  FORMULATION 


Formulation  of  the  delamination  problem  is  based  on  the  theory 
of  anisotropic  elasticity  for  nonhomogeneous  solids.  Well  known 
Lekhnitskii's  stress  functions  [13]  are  introduced  to  establish 
governing  partial  differential  equations  for  field  variables.  An 
eigenfunction  expansion  method  is  employed  for  determining  the  stress 
singularity  at  the  delamination  crack  tip.  The  boundary  collocation 
technique  is  then  used  to  evaluate  the  complete  solution  for  finite¬ 
dimensional  composite  laminates. 

2.1  Assumptions  and  Delamination  Model 

Consider  a  composite  laminate  (Fig.  1)  composed  of  unidirec¬ 
tional  fiber-reinforced  plies  of  uniform  thicknesses,  h^,  . h^. 

For  simplicity  but  without  loss  of  generality,  we  restrict  ourself  to 
the  cases  of  symmetric,  angle-ply  composite  laminates  with  fiber 
orientations  (Q . /Q ^) .  Ply  thicknesses  are  also  assumed  to 
be  symmetric  with  respect  to  the  x-z  plane,  i.e.,  h^  =  h^,  h^  = 

....  The  composite  has  a  finite  width  2b  and  is  subjected  to  a  uniform 
axial  extension,  e,  where  s  =  constant,  along  the  z-axis.  The  composite 
laminate  is  sufficiently  long  that,  in  the  region  far  away  from  the 
ends,  end  effects  are  negligible  by  virtue  of  the  Saint  Venant 
principle.  Consequently,  stresses  in  the  laminate  are  independent  of 
the  z-axis.  The  case  in  which  stresses  and  displacements  are  independ¬ 
ent  of  z  corresponds  to  the  well  known  generalized  plane  deformation 
[13].  Edge  delamination  occurs  in  the  form  of  a  crack  along  the 
interface  of  dissimilar  plies  with  fiber  orientations,  9^^  and 
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Perfect  bonding  is  assumed  in  the  composite  everywhere  except  the 
region  of  delamination. 

2.2  Basic  Equations 

For  each  individual  lamina,  the  constitutive  equations  in  the 
structural  coordinates  x-y-z  may  be  expressed  by  the  generalized  Hooke’s 
law  in  the  contracted  notation  as 


£.  =  S. ,  a. 
1  ij  1 


(i»j  ~  1,2,....6) 


(1) 


where  e.  and  a.  are  strain  and  stress  tensors,  and  S..,  the  compliance 
1  J  ij  ’ 

tensors,  respectively.  The  strain-displacement  relationships  are  given 
by 


=  e  =  u,  , 
1  X  ’x 


=  £  =  w,  , 

3  z  z 


^2  “  ^y  '  "^’y’ 


£  =  Y  =  +  V,  , 

4  '  yz  T  ’  z 


£c  =  Y  =  w,  +  u,  ,  e^  =  Y  =  u,  +  V,  , 
5  xz  X  z  6  xy  y  x 


(2a-b) 


(2c-d) 


(2e-f) 


where  the  subscript  after  comma  denotes  the  partial  differentiation  with 
respect  to  the  variable.  Under  these  assumptions,  there  remain  three 
comparability  relations : 


=  Y 


X 


yy 


XX 


xy; 


xy 


^“'Yxz,y  ^yz,x^’y 


(-y  +  Yxz  =  0- 


(3a) 

(3b) 

(3c) 


For  a  symmetrical  angle-ply  composite  laminate  subjected  to  uniform  axial 
strain  e,  it  can  be  shown  easily  that 
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(4) 


since  the  relative  artgle  of  rotation  for  the  symmetric  composite  lamin¬ 
ate  about  z-axis  vanishes . 

Based  on  the  definition  of  the  problem,  i.e.,  e  =  e  =  constant, 

z 

in  Eq  1  can  be  expressed  in  terms  of  other  stress  components  by 

'^z  "  "  ^3j  =  1,2, 4, 5, 6).  (5) 


Thus,  the  generalized  Hooke’s  law  may  be  modified  to  have  the  following 
form: 

(i,j  .  1.2,4.5,6),  (6a) 

where 


S.  , 


^3iSj^^33’ 


e . 

1 


i3'  33 


(i.j  ^  3). 


(6b) 


It  can  be  seen  from  Eq  6  that  e^  has  the  role  of  initial  strains  in  the 

laminate.  For  the  current  stress  formulation,  it  may  be  more  convenient 

to  introduce  the  initial  stress  a.  such  that 

JO 

4l  =  Sy(a.  -  (I.J  -  1.2, 4, 5, 6),  (7) 


where  a .  =  -S . .  e . . 

lo  ij  j 

into  two  parts,  i.e.. 


It  is  possible  to  deccm.pose  the  complete  solutions 


a . 
1 


a'"' 

1  J 


(8a) 


e . 

1 


where 


and  =  S..Co^P^-  a.  ). 

1  iJ  J  Jo" 


(8b) 


(9a-b) 
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3.  GENERAL  SOLUTIONS  AND  ASSOCIATED  STRESS  SINGULARITIES  FOR  DELAMINATION 


3.1  Solutions  for  and  uf^^ 

_ 1 _ 1 

Introducing  the- stress  functions  F(x,y)  and  '?'(x,y)  which  satisfy 
equations  of  equilibrium  identically  and  following  the  procedure  by 
Lekhnitskii  [13],  we  obtain  a  pair  of  coupled  partial  differential 
equations  as  follows: 


L^F  +  L^^  =  0, 


L^F  +  =  0, 


(10a) 


(10b) 


where  L^,  L^  and  L^  are  linear  differential  operators  of  the  second, 
third  and  fourth  order,  respectively,  defined  by 


T  =  5  3^  .  o  3^ 

2  44  ^  "  45  3x3y  55  3^  ’ 


(10c) 


S  '^24  ^^25  ^46^  3x'^3y  "  ^^14  ^56^  3x3y^ 


15  V  ’ 


(lOd) 


H  ^22  3x^  "  ^^26  3x3 3y  ^^^12  ^66^  3x^¥^ 


^^16  3x3y^  '''  ^11  37^  * 


(lOe) 


Lekhnitskii  [13]  has  shown  that  the  general  solution  for  Eqs  lOa-lOb  may 
be  expressed  as 


F(V  =  I  \(V’  ■ -I  Vi<V' 


(lla-b) 


where  Z,  =  x  +  y,y;  the  prime  (')  denotes  differentiation  of  the  function 
with  respect  to  its  argument;  and  are  the  roots  of  the  algebraic 
characteristic  equation 


i^(y)£2(y)  -  =  0,  (12a) 

and 

with 


■  ®44’ 

(13a) 

«3(P)  =  S^5»’ 

-  (Su-"  556)f"  +  (=25*  h6'>'‘  -  "24' 

(13b) 

-  “is'-’  ^  (2S33+  -  2S3^p  + 

(13c) 

Introducing  the  following  form  for  the  function 


5+2 

F^(Z^)  =  C^Z^  /[(5  +  2)(6  +  1)]  , 


(14) 


where  C,  and  5  are  arbitrary  complex  constants  to  be  determined  later, 

lx 


we  can  obtain  the  stress  and  displacement  expressions  as  follows ; 

=  y  rc  u^Z'^  +  C  • 

^x  ,  S  ^  Vk\  ^k+3  '^k  ’ 

k=l 


(ISa) 


(h)  _  V  Tf.  7'5  ,  p  ys, 

y  ■  k=i  ^ V  ’ 


(15b) 


vL  ^^k’^k^  \+3  ^k  ’ 

K-^X 


(15c) 


~  ^Vk’^k^  ^k+3  ‘^k  ’^k  V’ 

k=l 


(15d) 
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(15e) 


=  -I  [C  U  zj  +  C  u  Z?^’ 
xy  k  k  k  k+3  k  k 


and 


=  J  'VkC  "k+3  '>k  1> 


,5+1 


^+l3 


.(^i)  _  Y  rr  n  7^'^^  4.  r  TT  4.  1^ 

k=l 


^  fVkC'  ^V4-.  ?"■'■]/ (5  +  1) 

k=l 


,5+1  .  ^  —  -^+1-, 

"k+3  k  k 


where 


Pk  ~  ^ll’^k  ^12  ■  ^14\  ^IsVk  ■  ^16^k 


^k  *  ^12'^k  ^22'^^k  “  ^24'^k^^k  ■*'  ^25’^k  "  ^26  ’ 


\  “  ^14’^k  ■*■  ^24^^k  ■  ^44’^k^^k  ■'■  ®45^k  ■  ^46  * 


(16a) 

(16b) 

(16c) 

a6d) 

(16e) 

(16f) 


The  constant,  5,  in  Eqs  15  and  16  may  be  chosen  that  the  stresses 
and  displacements  and  satisfy  interface  continuity  and 

homogeneous  boundary  conditions.  Taking  complex  conjugate  of  Eqs  15 
and  16,  their  forms  are  invariant;  thus,  5  appears  as  a  set  of  complex 
conjugates,  which  enables  to  make  Eqs  15  and  16  real  functions  by  super¬ 
position.  Furthermore,  finiteness  of  displacements  at  the  origin  re¬ 
requires  that  Re[5]  >  -1,  where  Re  represents  the  real  part  of  6. 


3.2  Solutions  for  and  uf^^ 


.(p) 


Since  a,  and  e.  in  Eqs  6  and  7  are  constant,  we  may  choose  a 
10  1  ^  1 

in  Eq  9  as  constants  so  that  they  satisfy  the  equations  of  equilibrium 
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and  compatibility  conditions  identically.  For  each  individual 
let  take  the  following  form: 


(P)  _ 


X 


k=l 


(P)  _ 


+  y  (d,  +  d,  ) , 

yo  k  k^’ 


,(P)  _ 


yz 


““k 


,(P) 


xz 


^kZO  Ji'Vkl'k  +  ““k  \  >'k>> 


r<P>  - 

xy 


^xyo  ^  ‘^k  '^k^  * 


lamina, 

(17a) 

(17b) 

(17c) 

(17d) 

(17e) 


Substiuting  Eqs  17a-17e  into  Eq  9b  and  integrating  the  strain-displace¬ 
ment  relations,  we  obtain 


(  )  ^ 

*  ‘‘k  \  +  "3*  ■  “l"  +  ''o’ 

(  )  ^  _  _  _ 

w  ^  =  I  ^\Pk\  '*'  \  ^k  ~  “2^  ^  ‘*^1^  ^o’ 

k— 1 


where  u^,  v^,  w^  and  are  related  to  rigid-body  displacements  and  rota¬ 
tions.  The  complex  constants  d^^  are  required  to  satisfy  the  near-field 
traction  boundary  conditions  and  continuity  conditions  along  the  ply 
interface. 
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3.3  Delamination  Crack-Tip  Stress  Singularity 

Consider  a  delamination  between  two  plies,  say  the  kth  and  (k+l)th 
ply,  in  a  composite  subjected  to  general  loading  as  shown  in  Fig.  2. 
Assuming  that  interl^inar  crack  surfaces  are  free  from  traction,  we 

introduce  the  following  boundary  conditions  for  the  eigen  stresses 

(h) 

a. 

1 

=  0  (i  =  k,  <))  =  it;  i  =  k+1,  <(>  =  -tr).  (19) 

y  yz  xy 

The  superscript  h  is  dropped  in  the  expression  for  convenience. 

Continuity  conditions  for  displacements  and  interlaminar  stresses 
along  the  interface,  <f>  =  0, 


y  yz  xy  y  yz  xy 

Substituting  Eqs  15  and  16  into  19  and  20,  we  obtain  the  following 

(k)  (k+1) 

twelve  linear  algebraic  equations  in  C^  and  C^  : 


(20a) 


(20b) 


y  -  0, 

^  m  jm  m+i  jm 

m=l 


(j=l,2,3) 


(21a) 


y  =  0  (j=l,2.3), 

m  im  m+3  jm 

Q=1 

I  ._(k)„(k)  ,  „(k)  ^(k),  _  I  .  (k+1)  (k+1)  (k+1)  jr(k+l)^ 

^/^m  ^rm  ^  \+3  m  ^  ”  -/m  m  ^  V+3  m 


(21b) 


(r=l,2,3,4,5,6). 


(21c) 


,  „(k)_,  Ck)_  (k)  (k)_  (k)  j,(k)^  (k)  pCk)^  (k)  ^ 

where  =1,  p^  ,  q^ 

p(k)^  j^(k)^  Solving  Eqs  21c  for  C^^^  and  substituting  the  resulting 
6m  m 


4. 


expressions  into  Eqs  21a,  we  get 


I  I  I 


mn  jm  (m+3)n  jm  ^ 


(21d) 


Equations  21b  and  21d  consist  of  a  system  of  six  linear  homogeneous 

algebraic  equations.  For  nontrival  solutions  of  the  determinant 

m 

of  coefficients  of  the  algebraic  equations  must  vanish.  This  leads  to  a 
characteristic  equation  of  the  following  form: 

1)3[a(6)|  =  0,  (22) 

where  [a(6)|  is  a  3  by  3  determinant  involving  5  in  a  transcendental  form 

and  material  constants,  u  and  of  the  adjacent 

mm  m  m 

layers.  Details  of  A (5)  may  be  found  in  Ref.  [14].  The  general  form  of 
6,  which  are  the  eigenvalues  of  the  problem.,  may  be  written  as 


5n  =  n,  or  5^  =  (n  -  j)  ±  iy 


(n=0 ,1,2,...), 


where  y  is  a  constant  related  to  elastic  constants  of  adjacent  plies. 
Thus,  for  each  we  have  the  eigenfunctions  of  the  form  Eqs  15  and  16 
whose  coefficients  may  be  determined  from  the  remote  boundary  conditions 
other  than  Eq  19.  It  is  important  to  note  that  the  5^  bounded  by 


0  >  Re[5^]  >  -1 


characterize  the  inherent  stress  singularities  of  the  delamination  crack 
stresses  in  a  composite  laminate. 

For  cross-ply  composite  laminates,  the  differential  operator 
vanishes  identically.  Thus,  F(Z)  and  'l'(Z)  are  uncoupled,  and  the  form 
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of  Eq  23  can  be  simplified  and  expressed  explicitly  as 

5  =  n,  or 

n 

(n  -  ±  ^  £n{[b  +  (b2  -  4a2)^]/(2a)},  (25) 

~  (k)  'fk+l')  ('k') 

where  a  and  b  are  related  to  material  constants  S;.,  S;.  ,  u.  , 

ij  xj  ’  1  ’ 

-(k)  ,  (k+1)  -(k+1)  , 

and  ,  y^  shown  in  Appendix  1.  In  a  limiting  case  of 

isotropic  materials,  it  can  be  shown  that  6  have  the  form, 

n 

5^  =  Cn  -  ±  -1^  -  4v^^^)] 

(26) 

where  G  and  v  denote  the  shear  modulus  and  Poisson’s  ratio,  respectively. 
The  eigenvalues  of  Eq  26  were  first  obtained  by  William  [1]  and  later  by 
Zak,  at  al.  [2]  for  interfacial  cracks  in  isotropic  media. 
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4.  DELAMINATION  STRESS  INTENSITY  FACTORS  AND  ENERGY  RELEASE  RATES 


The  eigenfunctions  and  the  unknown  constants  for  and  in 

1  i 

Eqs  15  and  17  are  determined  by  imposing  appropriate  (materials  and 

geometric)  symmetry  and  traction  boundary  conditions,  which  will  be 

discussed  later.  Hence,  complete  stresses  and  displacements  and 

(a )  . 

u^  in  the  a-th  lamina  can  be  fully  established.  Neglecting  the 
higher-order  terms,  we  note  that  the  typical  structure  of  near-field 
crack  tip  stresses  can  be  shown  to  have  the  following  form: 


n  3 


where  ®ijk  related  to  the  material  constants,  geometry, 

and  boundary  conditions;  6^  are  eigenvalues  bounded  by 

-1  <  Re[5j]  <  0  to  insure  the  positive  definiteness  of  strain  energy 

of  the  elastic  body.  It  is  clear  that  the  eigenvalues  which  satisfy 

Eq  24  lead  to  as3rmptotic  near  field  stresses.  For  the  convenience  of 

(a) 

further  development,  the  stress  around  the  crack  tip  may  be  re¬ 

written  as 

(ct )  ^  (ct ) 

'^1-  ~  1  ^  n  A  (xjy;5.)  +  0  (non-singular,  higher-order  terms), 

j=l  ^  ^ 


where  is  the  j-th  singular  component  of  the  stress  cor¬ 

responding  to  the  eigenvalue  6^  which  meets  Eq.  24. 

In  the  context  of  mechanics  of  fracture,  it  is  possible  to  define 
the  so-called  stress  intensity  factors  for  the  delamination  in  a  manner 
analogous  to  that  given  in  Refs.  [4 >6]  by  considering  the  interlaminar 


stresses  ahead  of  the  crack  tip  along  the  interface,  i.e.. 
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(29a) 


II  _ 

K  =  lim  ^  /ItT  X  (x,0;6.), 

^  xr^-0+  j=l  ^  ^ 


-5  • 

K_t.  =  lim  I  /ItT  X  2a,  (x,0;6  ), 
x*0+  j-l  ^ 


K  =  lim  I  'PH  X  2a  (x,0;6  ) 
HI  4j  j 


(29b) 

(29c) 


where  the  superscript  a  is  omitted,  because  tractions,  and  a^, 

are  continuous  across  the  interface. 

While  the  stress  intensity  factors  and  describe  the 

details  of  the  delamination  crack-tip  field,  the  strain  energy  release 
rate  G  is  also  of  significant  interest,  since  this  is  a  quantity 
physically  measurable  in  experiments  and  mathematically  well  defined. 
The  fracture  energy  release  rate  in  a  delaminated  composite  may  be 
evaluated  by  using  Irwin’s  virtual  crack  extension  expression  [15], 


G  = 


'III 


6  3 

=  lim  ^  f  {a  (r,0)[v^^^6B-r,TT)  -  (6B-r  ,-Tr)  ] 

6S^0  Jo  y 

+  T^y(r,0)[u^^\5B-r,ir)  -  u^^‘^^\sB-r ,-Tr) ] 

+  T  (r  ,0)  [w^^^  C6B-r  ,u)  -  (66-r,-TT)]}dr  »  (30) 

yz 


where  polar  coordinates  (r,<)))  are  used  for  the  convenience  of  computa¬ 
tion.  The  interlaminar  stresses,  o^,  1-9  I*® 

obtained  from  the  crack-tip  stress  field  equations  such  as  Eq  27.  The 
corresponding  displacements  are  also  those  of  the  crack-tip  field 
equations  obtained  in  the  previous  section. 
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5.  NUMERICAL  EXAMPLES  AND  DISCUSSION 


The  formulation  and  analysis  for  the  problem  outlined  in  previous 
sections  have  been  programmed  into  a  solution  scheme  suitable  for 
numerical  computation.  For  the  purpose  of  illustrating  the  fundamental 
behavior  of  the  delamination  fracture  in  composite  laminates,  graphite- 
epoxy  systems  with  symmetric  C0/-0/-0/9)  fiber  orientation  containing 
edge  delamination  cracks  along  the  0  and  -0  ply  interface  are  studied. 
The  particular  material  system  and  ply  orientations  are  selected  here 
because  they  have  been  previously  Investigated  in  some  detail. 

The  composite  laminate  is  subjected  to  a  uniform  axial  extension  and 
has  a  geometry  shown  in  Fig.  1  with  a  width-to-thickness  ratio  2b/2W 
and  uniform  ply  thickness  h^.  Delamination  cracks  of  length  cc  are 
assumed  to  emanate  from  the  edges  of  the  composite.  Lamina  properties 
typical  of  high-modulus  unidirectional  graphite-epoxy  composite  for 
aircraft  construction  are  used  in  the  computation  (Table  1) .  For  com¬ 
posite  laminates  with  the  aforementioned  laminate  geometry  and  ply 
orientations,  several  geometric  and  material  symmetry  conditions  may 
be  introduced  to  simplify  the  formulation  further.  The  problem,  there¬ 
fore,  can  be  solved  very  conveniently  and  accurately. 

5.1  Symmetry  and  Boundary  Conditions,  and  Further  Simplifications 
The  symmetric  ply  orientations  and  geometry  of  the  composite 
laminate  (Fig.  3)  lead  to  the  following  conditions  for  displacements: 
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on  X  =  b-a. 


(31a) 


3u  _  Iw  _  ^ 

3y  9y  3x 

9v  3w  3u  -  , 

—  =  ^  =  -^=0  ony  =  -h_,  (31b) 

3x  3y  3y  ■'2’ 

where  the  origin  of  .the  coordinates  is  moved  to  the  left  tip  of  the 

delamination.  The  traction-free  boundary  conditions  on  edges  and  lateral 
surfaces  of  the  composite  laminate  may  be  written  as 


a=T  =T  =0  onx=  -a,  (31c) 

X  xy  xz 


a  =  T  =  T  =  0 
y  xy  yz 


on  y  =  h^. 


(31d) 


Thus,  only  a  quarter  of  the  laminate  cross  section  needs  to  be  considered. 
The  boundary  conditions,  Eqs  31a-d,  contain  arbitrariness  of  rigid  body 
displacements,  which  is  a  characteristic  of  traction  boundary  value 
problems . 

(h  ^ 

Since  the  eigen  solutions,  ,  satisfy  Eq  19  and  interface  con¬ 
tinuity  conditions  Eq  20,  we  require  that  satisfy  these  conditions 

(o.  ^ 

also.  To  determine  d^  uniquely  for  Eqs  17  and  18,  we  further  require 
the  particular  solutions  satisfy  the  following  conditions: 


u(P> 

=  0 

on 

X  =  b-a. 

(32a) 

v^P) 

=  0 

on 

y  =  -h2, 

(32b) 

w^P^ 

=  0 

at 

(0,0,0). 

(32c) 

Substituting  Eqs  17  and  18  into  Eqs  19,  20  and  32  gives 

(ci=l,2)  (33a) 

k=l 
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=  I  [d 


k  '^k 


(33b) 


»<“>  =  0, 
O 


(33c) 


(a)  (a) 


(33d) 


k=l 


(33e) 


I  [d,f'  +  d/“>l  .  -a<“> 


(34a) 


^2^[dk  +  d^ 


(34b) 


I  L\  +  dv  Wv  ]  = 


k  ^k 


(34c) 


I  [d(»pa).df>7<i>i  =  [  (d‘^p.dP;«b. 


k  "^k 


'k 


k  ‘^k 


(34d) 


I  [d(«t(».d«>T<ib  .  I  [d«)t«)  .  d«)  t<^>, , 


k  k 


k  k 


k  k 


(34e) 


V  r4(°‘)  4.  4.7^“)<'7(“)  4.  7(“)  7(“)m  -  n  /'■5/^^ 

L  [\  +  Pv  Pt  )  +  d^  (q^  +  )]  =  0.  (34f) 


■k  ^k 


k  ^^k 


k  '^k 


Equations  34a-f  give  ten  linear  algebraic  equations  for  twelve  real 
unknowns  for  3-nd  uf^\  Hence,  we  may  set 


Imtd^^^]  =  0 


(35a) 


to  reduce  the  additional  degrees  of  freedom.  For  symmetric  angle-ply 


laminates,  it  can  be  shown  that  Eq  34f  is  satisfied  identically. 
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Therefore,  instead  of  using  Eq  35a,  it  maybe  required 


(35b) 


Since  the  complete  solutions  for  stresses  and  displacements  must 

satisfy  the  symmetry  and  remote  boundary  conditions,  Eqs  31a-d ,  the 

(h) 

following  relations  can  be  established  immediately  to  evaluate  and 

„<w. 


+  afP)  .  0 
1  1 


(i=l,5,6)  on  X  =  -cz, 


(36a) 


„(h)  +  ^(P)  .  0 


(i=2,4,6)  on  y  =  h^. 


(36b) 


- =°’3r~  +  ai  X  =  b-a. 


(36c) 


— =°’ar“'^a^^ — =°’a7"+a?  ’  ^  = -^2* 


(36d) 


By  using  the  eigenfunctions  derived  previously  and  the  boundary  collocation 
method,  the  boundary  conditions  given  in  Eqs  36a-d  can  be  matched  con¬ 
veniently  in  the  least- square  sense.  Thus,  the  eigen  solutions  for 

0 and  uf^^  can  be  determined.  Numerical  solutions  for  the  problem 
1  1 

by  using  the  collocation  procedure  are  related  to  the  truncation  of 
eigenfunctions  and  number  of  collocation  stations.  Due  to  space 
limitation,  the  detailed  discussion  of  solution  convergence  and  accura¬ 
cy  is  reported  elsehwere  [16]  •  The  results  presented  in  this  section 
are  from  collocation  calculation,  which  has  a  maximum  mismatch  within 
one  percent  deviation  from  prescribed  boundary  conditions. 
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5.2  Stress  Singularity  for  Delamination  in  Composites 

Now  consider  a  delamination  lying  between  9  and  -6  plies  (Fig.  1) 
in  a  graphite-epoxy  composite  with  ply  properties  given  in  Table  1.  The 
interface  continuity  and  traction  boundary  conditions  along  crack  sur¬ 
faces  lead  to  a  standard  eigenvalue  problem  for  the  homogeneous  solution, 
as  discussed  in  Section  3.3.  The  eigenvalues  5^  obtained  from  the 
transcendental  equation  provide  basic  structures  of  near-field  stress 
and  displacement  solutions  for  the  delamination  problem.  The  order  of 
stress  singularity  and  the  asymptotic  nature  of  the  crack  tip  stresses 
depend  on  the  values  of  5^,  which  satisfy  the  constraint  condition  of 
Eq  24.  Thus,  the  eigenvalues  corresponding  to  this  restriction  are  of 
fundamental  importance  in  understanding  the  delamination  failure  be¬ 
havior.  For  edge  delaminated  (±9)  graphite-epoxy  composites  with 

s 

various  fiber  orientations,  the  eigenvalues  6^,  which  satisfy  the 
aforementioned  constraint  condition,  are  found  by  the  present  eigen 
analysis  and  given  in  Table  2.  The  stress  singularities  for  an  inter¬ 
face  crack  between  two  highly  anisotropic  laminae  are  observed  to  con¬ 
tain  a  pair  of  complex  conjugates,  2  ~  ~0*5±iY,  and  a  constant, 

5^  =  -0.5.  This  situation  is  unique  and  different  from  that  of  an 
interface  crack  between  two  isotropic  or  orthotropic  media  in  the  sense 
that  5^,  5 2  and  5^  exist  simultaneously  in  the  present  delamination 
problem.  In  the  degenerated  cases  such  as  ±9  =  0°  and  90°,  the  composite 
laminates  become  unidirectional.  The  delamination  is  located  in  an 
ortho tropic  material;  the  classical  inverse  square-root  singularity 
for  crack-tip  stresses  is  recovered  fully.  It  is  noted  that  the 
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present  physical  model  and  the  eigenfunction  analysis  lead  to  an  oscilla¬ 
tory  stress  singularity,  as  are  the  cases  of  interface  cracks  in  iso¬ 
tropic  or  orthotropic  materials. 


5.3  Asymptotic  Stress  Field  Around  Delamination 

Complete  solutions  for  delamination  cracks  in  finite  dimensional 
composite  laminates  are  obtainable  by  using  the  present  Lekhnitskii’s 
complex  stress  potential  formulation  and  eigenfunction  expansion.  With 
the  aid  of  the  boundary  collocation  method,  the  asymptotic  stress  field 
around  a  delamination  may  be  expressed  in  a  general  form  as 


a . 

1 


k=l 


jk 


-0 . 5+iY 
'k 


+  D 


j (k+3) 


_-0 . 5+iY 
Z  )  + 


-0.5-iY 

'k 


(k+3) 


_-0.5-iY 
Z  ) 


+  (F 


,-0.5 


jk  ‘-k 


(k+3) 


_-0.5 
Z  )} 


(j=l»2, . ,6) , 


(37) 


where  D.,  ,  E.,  ,  and  F.,  are  known  quantities  satisfying  the  following 
j  k  j  k  j  fc 

relations : 


°jk  ^j(k+3)’  °j(k+3)  ^jk’ 


^jk  ^j(k+3)’ 


(38) 


to  insure  being  real.  More  concisely,  cr^  can  be  written  as 


a.  =  r  2  [A.cos(y  r)  +  B.sin(Y  r)  +  C.]. 
3  J  3  3 


(39) 


For  illustrative  purposes,  the  structures  of  near-field  stresses  and 
displacements  ahead  of  a  delamination  (r,0)  are  given  for  a 
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(45°/-45°/-45°/45'’)  graphite-epoxy  laminate  with  h^  =  h2  =  1  in., 
cz  =  1  in.  and  2b/h  =  4  subjected  to  =  e  as  follows; 

c?y(r,0)  =  [0.04339cos(0. 03434  S.n  r)+0.39498sin(0. 03434  In  r)] 

+  0(1),  (40a)  ' 

-0.5 

T  (r,0)  =  [0.45347  cos(0. 03434  iln  r)-0.04981sin(0. 03434  iin  r)]  r 
yz 

+  0(1),  (40b) 

T^y(r,0)  =-0.002449  r  +  0(1)  ,  (40c) 

and 

(6P-r,Tr)  =  {-0.29576cos[0. 03434  £n(56-r)] 

+  0.01208sin[0. 03434  Jln(68-r)]}  (5B-r)°'^ 

-  0.003403(d6-r)°*^  +  0(1),  (41a) 

(5B-r,TT)  =  {0.02888cos[0. 03434  Jln(63-r)] 

+  0.70682sin[0. 03434  Jin(53-r)]  }  (63-r)°'^ 

+  0(1),  (41b) 

w^^^(S3-r,u)  =  {0.61564  cos[0. 03434  !ln(5B-r)] 

-  0.025155  sin[0. 03434  £n(5B-r) ] } (66-r)° 

+  0.001602(66-r)°*^  +  0(1),  (41c) 

where  the  components  of  stress  are  scaled  by  10®  x  e^,  and  the 
displacements  by  It  is  noted  that  the  elastic  stresses  near  the 
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delamination  crack  tip  in  a  composite  laminate  possesses  the  well  known 
oscillatory  behavior,  and  the  displacement  field  also  exhibits  an 
oscillatory  nature  with  crack  surfaces  overlapping  each  other.  As 
first  pointed  out  by  Malyshev,  et  al.  [7]  and  later  by  England  [5]  and 
Erdogan  [8]  for  interfacial  cracks  between  dissimilar  isotropic  media, 
the  phenomenon  of  crack  surface  overlapping  is  confined  to  an  extremely 
small  region  and  the  interpenetration  is  not  of  significance  in  prac¬ 
tical  terms  of  fracture  mechanics.  However,  for  certain  combinations 
of  material  properties,  ply  orientations  and  loading  conditions  in 
composite  laminates,  the  crack  surface  contact  region  has  been  found  to  be 
extremely  large  [16]  .  Thus  the  current  model  needs  to  be  modified  to 
account  for  the  crack  surface  closure  and  contact  stresses  [16]  .  Studies 
on  interface  crack  closure  in  dissimilar  isotropic  media  were  reported 
recently  by  Comninou  [17],  Atkinson  [18]  and  Achenback  [11]. 

5.4  Delamination  Crack  Tip  Stress  Intensity  Factors 

Since  the  Irwin  fracture  criterion  is  local  in  nature  and  requires 
precise  knowledge  of  the  local  conditions  at  the  delamination  crack  tip, 
the  stress  intensity  solutions  are  obviously  of  great  significance. 
According  to  the  present  fracture  mechanics  theory  of  composite  delami¬ 
nation,  stress  intensity  factors,  K^,  and  ™3y  be  evaluated  by 

the  rigorous  analysis  described  in  Section  4.  The  K^,  and 
lead  to  detailed  information  of  the  stress  and  displacement  fields  in 
the  neighborhood  of  the  delamination  crack  tip,  and  may  relate  to  the 
onset  of  delamination  extension  upon  reaching  a  critical  level.  The 
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magnitudes  of  shown  in  the  formulation  depend  on  the  delamination 
length,  ply  orientations,  laminate  geometry,  and  loading  conditions. 

Consider  the  (9/-0/-0/0)  graphite-epoxy  composites  with  various 

fiber  orientations  subjected  to  uniform  axial  extension  e  •  For  illus- 

z 

trative  purposes,  we  choose  a  composite  with  a  width-to -thickness  ratio 
2b/ (2h)  equal  to  8,  ply  thickness  h^  =  h^  =  1  in.,  and  delamination 
length  a  =  1  in.  The  mixed  mode  K^,  and  3.V&  determined  and 

given  in  Table  3  for  various  0*s.  It  is  observed  that  even  though  the 
composite  laminate  is  under  the  simplest  loading  condition,  the  delamina¬ 
tion  crack  tip  response  is  very  complicated  due  to  the  complex  inter¬ 
laminar  stress  distribution,  the  nonhomogeneity  of  the  solid,  the  aniso¬ 
tropic  ply  properties,  and  the  unusual  delamination  configuration  with 
respect  to  the  loading  direction.  The  out-of -plane  tearing  mode  stress 
intensity  factor  caused  by  interlaminar  shear  is  about  one  or 

two  orders  of  magnitude  higher  than  and  in  general  in  the  lami¬ 
nates  studied.  The  opening  mode  stress  intensity  is  also  very  sig¬ 
nificant  due  to  the  interlaminar  normal  stresses  cr^.  The  simultaneous 
presence  of  K^,  and  in  the  delamination  problem  is  unique  to 

angle-ply  fiber  composites,  and  is  not  observed  in  fracture  problems 
for  bonded  dissimilar  isotropic  media  in  general.  The  delamination  be¬ 
havior  is  inherently  three  dimensional  in  nature;  for  composites  with 
more  general  laminations,  crack  geometry  and  loading  conditions,  fully 
three-dimensional  stress  and  fracture  analyses  are  essential  for  ob¬ 
taining  complete  information. 

The  influence  of  laminate  geometric  variables  on  the  delamination 
behavior  is  best  illustrated  by  examining  the  changes  of  K^,  and 
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Kjjj  with  the  relative  thickness  of  upper  and  lower  plies  ^ 

(45°/-45‘’/-45'’/45'’)  graphite-epoxy  composite  (with  h^  +  h2  =  W  =  2  in.)* 
Given  the  crack  length,  (a  =  1  in.),  laminate  dimensions  (2b  =  4  in.), 
and  the  loading  condition  as  previously,  the  delamination  stress  inten¬ 
sities  for  various  h2^/h2's  are  shown  in  Fig.  4.  The  crack  tip  tear¬ 
ing  and  opening  stresses  and  have  a  maximum  intensification  as 
the  ply  thicknesses  h^  become  identical,  i.e.,  =  1.  The 

however,  reaches  a  minimum  due  to  the  reduction  of  t  .  It  should  be 
noted  here  that  the  K^,  and  depend  on  material  constants  of 

all  plies  as  well  as  the  overall  geometry.  Therefore,  the  dependency 
of  on  ply  properties  is  not  a  simple  matter  of  identifying  them  with 
geometric  variables,  and  they  may  not  have  the  simple  physical  inter¬ 
pretation  as  in  the  homogeneous  case. 


5.5  Strain  Energy  Release  Rates  for  Delamination 

The  equilibrium  and  stability  of  delamination  are  commonly  examined 
from  an  energy  rate  point  of  view.  The  strain  energy  release  rate  G  de¬ 
fined  in  Eq  30  is  a  quantity  characterizing  the  driving  force  for  delam¬ 
ination  extension.  The  delamination-growth  driving  force  can  be  easily 
determined  after  the  establishment  of  the  local  as3nnptotic  stress  and 
displacement  fields.  For  the  edge  delamination  problem  in  graphite-epoxy 
composites  considered  here,  the  G  value  may  be  obtained  in  a  general 
form  as 


G  = 


=  lira 
50^0 


1 

250 


{A  cos[y  - ^)]  +  A„sin[Y  £n(^-^ - ^)  ]  +  A  } 

1  r  z  r  J 


dr  , 


(42) 
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where  y  =  Ini[5^]  .  Equation  42  has  the  form  similar  to  the  one  derived 
previously  for  an  elastic-half  space  problem  by  Willis  [19] .  The  singu¬ 
lar  integration  may  be  carried  out  by  defining  an  anlytical  function 
with  the  cut  0  <  x  <  5S 


f(z)  =  (• 


z  -  63 


0.5+iY 


) 


(43a) 


so  that  we  have 

^^.63  -  rx^'^cos[Y  ^  ^)]dr  =  TT63/(e  +  e^^),  (43b) 

•  0 

f63  -o  0.5  to 

( - —)  sin[Y  ^'n( - ^)]dr  =  2  Y'n''S3/(e  +  e^^).  (43c) 

Jq  ^  r 


Thus,  the  total  energy  release  rate  G  can  be  determined  immediately  by 
substituting  Eq  43  into  Eq  42.  Table  4  shows  the  change  of  G  values  with 
ply  orientations  for  the  C0/-0/-9/e)  graphite-epoxy  composites  with  the 
material  properties  in  Table  1. 

To  study  the  basic  nature  of  delamination  extension  in  angle-ply 
composites,  strain  energy  release  rates  in  the  (45'’/-45'’/-45°/45“) 
graphite-epoxy  with  various  crack  lengths  are  examined.  Effects  of 
laminate  width  on  delamination  crack  extension  is  also  investigated. 

The  change  of  total  strain  energy  release  rate  G  with  delamination  length 
a  is  given  in  Fig.  5  to  illustrate  fundamental  characteristics  of  the 
delamination  fracture.  For  the  composite  laminates  with  various  2b/2h’s, 
the  G  is  observed  to  change  with  delamination  length  in  a  unique  manner . 
The  maximum  energy  release  rate  or  crack  extension  driving  force  occurs 
at  a  delamination  length  approximately  equal  to  one  or  two  ply  thick¬ 
nesses  in  the  composite  studied,  depending  on  the  (2b/2h)  ratio.  As 
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the  delamination  exceeds  this  characteristic  dimension,  G  decreases 
monotonically . 

On  the  basis  of  fracture  mechanics,  several  important  features 
regarding  delamination  fracture  are  revealed  from  the  Figure.  Assuming 
that  the  material  resistance  to  delamination  growth  remains  constant 
(i.e.,  the  failure  criterion,  =  constant,  is  used),  we  can  immedi¬ 
ately  conclude  that  there  exists  a  critical  delamination  length  asso¬ 
ciated  with  the  maximum  G  (for  example,  a*  -  2h  for  the  case  b/h  =  8) 
for  each  composite  laminate;  the  word  "critical"  means  the  one  that 
experiences  stable  crack  extension  at  the  lowest  load.  It  also  indi¬ 
cates  that  any  interlaminar  edge  flaw  inherently  in  the  composite, 
which  is  less  than  a  ,  will  experience  rapidly  unstable  growth  as  the 
load  or  G  reaches  a  critical  level,  and  is  anticipated  to  be  arrested 
at  a  later  stage.  Any  initial  delamination  greater  than  a'  will  ex¬ 
perience  a  stable  growth  under  monotonically  rising  loads;  that  is, 
there  exists  an  inherently  built-in  crack  arrest  mechanism  for  edge 
delamination.  These  phenomena  predicted  by  the  G  -  a  curve  have  been 
noted  by  several  researchers  conducting  experimental  and  analytical 
studies  on  the  delamination  fracture.  The  a  may  be  an  important 
quantity  in  the  life  prediction  for  delaminated  composite  materials 
and  structures  subjected  to  static  and  cyclic  loading. 

The  delamination  strain  energy  release  rate  is  also  a  function  of 
other  geometric  variables.  For  example,  G  is  significantly  affected 
by  the  relative  ply  thickness  h^/h2.  In  a  (45°/-45®/-45°/45‘’)  graphite- 
epoxy  with  a  geometry  given  before,  the  change  of  G  with  k^/h^  is  given 
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in  Fig.  6,  where  maximum  driving  force  occurs  at  indicating  the 

criticality  of  the  relative  ply  thickness  to  delamination  fracture  in 
composites . 

It  is  noted  here  that,  even  though  the  near-field  stresses  possess 
an  oscillatory  singularity  and  may  not  have  the  usual  significance 
attached  to  them  as  in  the  cohesive  (homogeneous)  case,  the  energy  re¬ 
lease  rate  G  is  well  defined  mathematically  and  physically,  and  should 
be  the  quantity  of  major  interest.  The  G  and  its  components  G^,  G^^  and 
Gjjj  can  be  evaluated  theoretically  and  experimentally  to  provide  a 
basic  measure  of  the  delamination  fracture. 
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6.  SUMMARY  AND  CONCLUSIONS 


An  analytical  method  for  studying  delamination  is  presented  in  this 
paper.  Fundamental  nature  of  edge  delamination  in  advanced  fiber  composite 
laminates  is  examined.  Based  on  the  theory  of  anisotropic  elasticity, 
the  composite  delamination  problem  is  formulated  by  using  Lekhnitskii's 
comp lex- variable  stress  potentials  and  an  eigenfunction  expansion  method. 
Exact  orders  of  the  three-dimensional  stress  singularity  in  a  delamina¬ 
tion  crack  tip  region  are  determined  from  the  eigen  analysis.  With  the 
aid  of  a  boundary  collocation  technique,  complete  stress  and  displacement 
fields  in  a  finite-dimensional,  delaminated  composite  are  fully  determined. 
Fracture  mechanics  parameters  such  as  the  mixed-mode  stress  intensity 
factors  and  associated  energy  release  rates  for  edge  delamination  are 
calculated  explicitly.  Solutions  are  obtained  for  edge-delaminated 
(9/-9/-0/9)  angle-ply  composites  under  uniform  axial  extension.  Effects 
of  delamination  lengths,  fiber  orientations,  lamination  and  geometric 
variables  are  studied  in  detail. 

Based  on  the  information  given  in  the  previous  sections,  the  follow¬ 
ing  conclusions  may  be  drawn: 

1.  An  analytical  method  based  on  the  theory  of  anisotropic 
elasticity  is  successfully  developed  to  study  edge  de¬ 
lamination  in  angle-ply  composite  laminates.  Formulation 
of  the  problem  is  carried  out  by  using  Lekhnitskii's  complex 
stress  functions.  Stress  singularities  for  delamination  between 
highly  anisotropic  laminae  are  obtainable  by  using  an  eigen¬ 
function  analysis .  The  order  of  delamination  crack-tip 
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stress  singularity  is  different  from  that  of  an  interface 
crack  between  dissimilar  isotropic  or  orthotropic  media 
by  the  simultaneous  presence  of  three  characteristic 
eigenvalues  of  -0.5+iY»  -0.5-iYj  and  -0.5. 

3.  The  fracture  mechanics  concept  may  be  extended  to  delamina¬ 

tion  problems  in  anisotropic  composite  laminates  by  properly 
defining  the  interlaminar  crack-tip  stress  intensity  factors 
such  as  Eqs  29a-c  and  strain  energy  release  rates.  For 
angle-ply  composite  laminates,  K^,  and  always 

occur  simultaneously  for  an  edge  delamination  with  being 

one  or  two  orders  of  magnitude  higher  than  the  other  two. 

4.  Complete  stress  and  displacement  fields  in  a  delaminated 
composite  may  be  accurately  determined  by  a  combined  eigen¬ 
function  expansion  and  a  boundary  collocation  method.  The 
asymptotic  solutions  are  characterized  by  (i=I,II,III) 
and  possess  the  well  known  oscillatory  behavior.  The  crack 
surface  overlapping  could  be  very  large  in  some  composite 
systems  with  certain  combinations  of  fiber  orientations,  ply 
stacking  sequences,  and  loading  conditions;  modifications 
[16]  of  the  current  model  to  include  crack  closure  may  be 
needed  for  these  cases . 

5.  The  crack  extension  driving  force  or  strain  energy  release 
rate  for  edge  delamination  in  composite  laminates  can  be  ac¬ 
curately  determined  by  using  Irwin’s  crack  extension 
concept.  Delamination  stability  in  composite  laminates 
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under  monotonically  rising  loads  can  be  assessed  for  any 


inherent  interlaminar  flaw  relative  to  the 
lamination  size  a*  obtained  in  the  current 


critical  de- 
G  -  a  curve. 
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TABLE  1 


Material's  Constarits  for  Graphite/Epoxy  Composite  Lamina 


20.0  X  10®  psi 
2.1  X  10®  psi 


^Lz  "  Sz 


V  =  V 

Lz  Tz 


0.85  X  10®  psi 

0.21 
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TABLE  2 


Dominant  Stress  Singularities  for  Delamination  in  (6/-0/-6/0) 

Graphite-Epoxy  Composites 


9 

^1.2 

CO 

0° 

-  0.5 

-  0.5 

15° 

-  0.5  ±  0. 006421 

-  0.5 

30° 

-  0.5  ±  0. 023991 

-  0.5 

45° 

-  0.5  ±  0. 034341 

-  0.5 

60° 

-  0.5  ±  0.029421 

-  0.5 

75° 

-  0.5  ±  0. 015791 

-  0.5 

o 

o 

-  0.5 

-  0.5 
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TABLE  3 


Stress  Intensity  Factors  K^*  for  Edge  Delamination  in  (9/-0/-G/0) 
Angle-Ply  Graphite-Epoxy  Composite"^  Subjected  to  Uniform  Axial  Strain 


±0 

fl 

V  tt 

15° 

0.08645 

0.01214 

-4.5588 

30° 

0.2330 

0.03330 

-3.6604 

45° 

0.1347 

0.01380 

-1.2968 

60° 

0.02025 

0.001360 

-0.1775 

75° 

0.006268 

-0.0002948 

0.0818 

K.  (psi  -  /in.)  are  scaled  by  10°  e 
1  z 

^a  =  bj^  ~  ^2  ~  ^  in.,  b  =  8  in. 

' 'For  the  delamination  crack  in  the  first  quadrant 
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TABLE  4 


Energy  Release  Rate  G  for  Edge  Delamination  in  (O/-9/-0/0)  Angle  -  Ply 
Graphite-Epoxy  Compositet  Under  Uniaxial  Extension 


±0 

G/10^ 

z 

15“ 

8.1076 

30“ 

4.0506 

45“ 

0.5740 

60“ 

0.0138 

75“ 

0.0036 

1  in. , 


b  =  8  in. 
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FREE  EDGE  INTERFACE 
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FIG.  1  EDGE  DELAMINATION  GEOMETRY  AND  LAMINATE  CONFIGURATION. 


FIG.  2 

EDGE  DELAMINATION  BETWEEN  kTH  AND  (k+l)TH  PLIES 
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4] 


SYMMETRY  AND  BOUNDARY  CONDITIONS  FOR  (e/-e/-e/0)  GRAPHITE-EPOXY  COMPOSITE 
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STRAIN  ENERGY  RELEASE  RATE  G  vs  DELAMINATION  CRACK  LENGTH  a  IN  (45°/-45°/-45°/45°)  GRAPHITE-EPOXY 
COMPOSITES  WITH  VARIOUS  LAMINATE  WIDTHS  h's  (hi  =  ho  =  1  In.) 


GRAPHITE /EPOXY  COMPOSITES  WITH 


APPENDIX  1 


Materials  Parameters  for  'Helamination  Crack-Tip  Eigenvalues 
in  Cross-Ply  Composite  Laminate 


5  =  n 

n 


or 


=  (n  -  ^)  ±  £n{[b  +  (b2  -  4a2)^]/(2a)} 

where 

a  =  -Mf -  Mf  *  Mf 

+ 

4  4 


and 


M, 


{“>  -  i  +  u^”))  -  (af)  +rf>)}. 


M, 


2“’  -2 


-(a)  ^(«)(„(«)  +„(«)),_ 


M, 


Cct)  _  ~  Ca)„  Ca)  Cot)  '  (a) 
3  ■  ^11  ^1  ^2  “  ^12 


M 


(a)  _~(a)~Ca)  _  ('c(“)')2  ,  m  7^“^') 

4  ^11^22  ‘■^12^  ^  ^11  ^^1  ’^2  ^^1  ’'2  ^ 


- 1  xf h(pf >  +  7’)] 
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